Linear dynamics of quantum-classical hybrids 
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A formulation of quantum-classical hybrid dynamics is presented, which concerns the direct cou- 
pling of classical and quantum mechanical degrees of freedom. It is of interest for applications in 
quantum mechanical approximation schemes and may be relevant for the foundations of quantum 
mechanics, in particular, when it comes to experiments exploring the quantum-classical border. The 
present linear theory differs from the nonlinear ensemble theory of Hall and Reginatto, but shares 
with it to fulfil all consistency requirements discussed in the literature, while earlier attempts failed 
in this respect. Our work is based on the representation of quantum mechanics in the framework 
of classical analytical mechanics by A. Heslot, showing that notions of states in phase space, ob- 
servables, Poisson brackets, and related canonical transformations can be naturally extended to 
quantum mechanics. This is suitably generalized for quantum-classical hybrids here. 
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I. INTRODUCTION 

The hypothetical direct coupling of quantum mechani- 
cal and classical degrees of freedom - "hybrid dynamics" 
- presents a departure from quantum mechanics that has 
been researched through decades for practical as well as 
theoretical reasons. In particular, the standard Copen- 
hagen interpretation has led to the unresolved measure- 
ment problem which, together with the fact that quan- 
tum mechanics needs such interpretation, in order to be 
operationally well defined, may indicate that it deserves 
amendments. In this context, it has been recognized 
early on that a theory which dynamically bridges the 
quantum-classical divide should have an impact on the 
measurement problem besides being essential for at- 
tempts to describe consistently the interaction between 
quantum matter and classical spacetime [2J. 

Numerous works have appeared, in order to formulate 
hybrid dynamics in a satisfactory way. However, they 
were generally found to be deficient for one or another 
reason. This has led to various no-go theorems accom- 
panying a list of desirable properties or consistency re- 
quirements, see, for example, Refs. 

• Conservation of energy. 

• Conservation and positivity of probability. 

• Separability of quantum and classical subsystems 
in the absence of their interaction, recovering the 
correct quantum and classical equations of motion. 

• Consistent definitions of states and observables; ex- 
istence of a Lie bracket structure on the algebra of 
observables that suitably generalizes Poisson and 
commutator brackets. 

• Existence of canonical transformations generated 
by the observables; invariance of the classical sector 
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under canonical transformations performed on the 
quantum sector only and vice versa. 

• Existence of generalized Ehrenfest relations {i.e. 
the correspondence limit) which, for bilinearly cou- 
pled classical and quantum oscillators, are to as- 
sume the form of the classical equations of motion 
( "Peres- Terno benchmark" test Q). 

These issues have been reviewed in recent works by 
Hall and Reginatto. Furthermore, there, they have in- 
troduced the first viable theory of hybrid dynamics that 
agrees with all points listed above [a-Q- Their ensemble 
theory is based on configuration space, which entails a 
certain nonlinearity of the action functional from which 
it is derived. This nonlinearity leads to effects and a 
poposal to possibly falsify the theory experimentally Q . 
We will comment on this issue in due course (Subsec- 
tion 5.3.). 

In fact, the aim of the present paper is to set up an 
alternative theory of hybrid dynamics, which is based 
on notions of phase space. This is partly motivated by 
work on related topics of general linear dynamics and 
classical path integrals [13, [Hi • Presently, we will extend 
the work of Heslot, who has demonstrated that quantum 
mechanics can entirely be rephrased in the language and 
formalism of classical analytical mechanics [12]. We thus 
introduce unified notions of states on phase space, ob- 
servables, canonical transformations, and a generalized 
quantum-classical Poisson bracket in particular. This 
will lead to an intrinsically linear hybrid theory, which 
fulfils all consistency requirements as well. 

It may be worth while to comment on the relevance 
of hybrid dynamics, even if one is not inclined to modify 
certain ingredients of quantum theory. There is clearly 
practical interest in various forms of hybrid dynamics, in 
particular in nuclear, atomic, or molecular physics. The 
Born-Oppenheimer approximation, for example, is based 
on a separation of interacting slow and fast degrees of 
freedom of a compound object. The former are treated 
as approximately classical while the latter as of quan- 
tum mechanical nature. Furthermore, mean field theory. 
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based on the expansion of quantum mechanical variables 
into a classical part plus quantum fluctuations, leads to 
another approximation scheme and another form of hy- 
brid dynamics. This has been reviewed more generally 
for macroscopic quantum phenomena in Ref. [13]. 

In all these cases hybrid dynamics is considered as an 
approximate description of an intrinsically quantum me- 
chanical object. Which can lead to new insights, for ex- 
ample, to an alternative derivation of geometric forces 
and Berry's phase fl4|. 

Such considerations are and will become increasingly 
important for the precise manipulation of quantum me- 
chanical objects by apparently and for all practical 
purposes classical means, especially in the mesoscopic 
regime. 

Furthermore, the backreaction effect of quantum fluc- 
tuations on classical degrees of freedom might be of con- 
siderable importance, in particular, if they originate in 
physically distinct ways. We recall continuing discussions 
of the "semiclassical" Einstein equation coupling the clas- 
sical metric of spacetime to the expectation value of the 
energy-momentum tensor of quantized matter fields. Can 
this be made into a consistent hybrid theory leaving grav- 
ity unquantized? This has recently been re-examined, for 
example, in Refs. d 1 [H, [3; various related aspects 
have been discussed, for example, in Refs. [T7l - [23| . 

Finally, concerning speculative ideas about the emer- 
gence of quantum mechanics from a coarse-grained deter- 
ministic dynamics (see, for example, Refs. [23 - l26j with 
numerous references to earlier work) the backreaction 
problem can be more provocatively stated as the problem 
of the interplay of fluctuations among underlying deter- 
ministic and emergent quantum mechanical degrees of 
freedom. Or, to put it differently: "Can quantum me- 
chanics be seeded?" 

The remainder of the paper is organized as follows. In 
Section 2., the results of Heslot's work are represented, in 
order to make the paper selfcontained, and we shall fre- 
quently refer to it in what follows. In Section 3., we intro- 
duce hybrid phase space ensembles and, in particular, the 
quantum-classical Poisson bracket which is central to our 
approach; the important issue of separability is resolved 
and time evolution discussed. In Section 4., hybrid dy- 
namics is studied, incorporating quantum-classical inter- 
action. Energy conservation, Ehrenfest relations, espe- 
cially for bilinearly coupled oscillators, are derived there. 
In Section 5., we discuss various aspects of the present 
theory, in particular, the possibility to have classical- 
environment induced decoherence, the quantum-classical 
backreaction, a deviation from the Hall-Reginatto pro- 
posal predicted by the hybrid dynamics developed in this 
paper, and the closure of the algebra of hybrid observ- 
ables. Section 6. presents concluding remarks. 



II. HAMILTONIAN DYNAMICS REVISITED 

In the following two subsections, we will briefly present 
some important results drawn from the remarkably clear 
exposition of classical Hamiltonian mechanics and its 
generalization incorporating quantum mechanics by Hes- 
lot [l2|. This will form the starting point of our discus- 
sion of the hypothetical direct coupling between quantum 
and classical degrees of freedom in Section 3. 

A. Classical mechanics 

The evolution of a classical object is described with 
respect to its 2n-dimensional phase space, which is iden- 
tified as its state space. A real- valued regular function on 
the state space defines an observable^ i.e., a differentiable 
function on this smooth manifold. 

Darboux's theorem shows that there always exist (lo- 
cal) systems of so-called canonical coordinates, commonly 
denoted by {xk,Pk), fc = 1, . . . ,n, such that the Poisson 
bracket of any pair of observables /, g assumes the stan- 
dard form 12711: 



/df_dg_ 



df dg \ 



'-'Pk 



(1) 



This is consistent with {xk,pi} = Ski, {xk,xi} = 
{Pk,Pi} = 0, k,l = l,...,n, and refiects the bilinear- 
ity, antisymmetry, derivation-like product formula, and 
Jacobi identity which define a Lie bracket operation, 
f,g — > {f,g}, mapping a pair of observables to an ob- 
servable. 

Compatibility with the Poisson bracket structure re- 
stricts general transformations G of the state space to 
so-called canonical transformations which do not change 
the physical properties of the object under study; e.g., 
a translation, a rotation, a change of inertial frame, or 
evolution in time. Such G induces a change of an ob- 
servable, / — >■ G{f), and is an automorphism of the state 
space compatible with its Poisson bracket structure, if 
and only if, for any pair of observables /, g: 



G{{f,9}) = {G{f),G{9)} 



(2) 



Due to the Lie group structure of the set of canonical 
transformations, it is sufficient to consider infinitesimal 
transformations generated by the elements of the corre- 
sponding Lie algebra. Then, an infinitesimal transforma- 
tion G is canonical, if and only if for any observable / the 
map f ^ G{f ) is given by / ^ /' = / -f^ {f,g}Sa, with 
some observable g, the so-called generator of G, and Sa 
an infinitesimal real parameter. 

Thus, for the canonical coordinates, in particular, an 
infinitesimal canonical transformation amounts to: 



Xk 



dg_ 
dpk 
dg 



5a 



Pk Pk ^Pk- TT—Sa , 

OXk 



(3) 



(4) 
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employing the Poisson bracket given in Eq. ([T]) . 

This analysis shows the fundamental relation between 
observables and generators of infinitesimal canonical 
transformations in classical Hamiltonian mechanics. 

For example, the energy Hcl given by the classical 
Hamiltonian function is the generator of time evolution: 
for g = Hch and 6a — St, the Eqs. ([3|) and (|4]) are equiva- 
lent to Hamilton's equations, considering an infinitesimal 
time step St. 



B. Quantum mechanics 



Finally, considering Hamilton's equations, deriving from 

n-. 

dm = a^=-*^l*) ' (9) 

m\ = -^ = *(n|i7 , (10) 

we see indeed that they represent Schrodinger's equation 
and its adjoint, using (H| — keeping the essential 

normalization condition ^ in mind. 



An important achievement of Heslot's work is the real- 
ization that the analysis summarized in the previous sub- 
section can be generalized and applied to quantum me- 
chanics; in particular, the dynamical aspects of quantum 
mechanics thus find a description in classical terms bor- 
rowed from Hamiltonian mechanics. Some related ideas 
have been presented earlier in Ref [28j . 



1. Preliminaries 

To begin with, we recall that the Schrodinger equation 
and its adjoint can be obtained by requiring the variation 
with respect to state vector \^) and adjoint state vector 
{^\, respectively, of the following action S to vanish: 



S:^ dt {^{t)\{idt - H)\^{t)) = dtL 



(5) 



which involves the self-adjoint Hamilton operator H per- 
taining to the physical object under study. - The adjoint 
equation follows after a partial integration, provided the 
surface terms do not contribute. This is guaranteed by 
the normalization condition: 



= constant = 1 



(6) 



which is an essential ingredient of the probability inter- 
pretation associated with state vectors. - Adding here 
that state vectors that differ by an unphysical constant 
phase are to be identified, we recover that the quantum 
mechanical state space is formed by the rays of the un- 
derlying Hilbert space, i.e., forming a complex projective 
space. 

Making use of the Lagrangian L, defined as the in- 
tegrand of the above action 5, we define a momentum 
conjugate to the state vector: 



(HI := 



dL 

d\¥) 



= im 



(7) 



with l^") :— dtl'i'), and obtain the corresponding Hamil- 
tonian function: 



(HI*) 



-i{n\H\^) =: ^(n,*) 



(8) 



2. The oscillator representation 

Quantum mechanical evolution can be described by a 
unitary transformation, |5'(t)) — U {t — to)\'^ {to)) , with 
U{t — to) = exp[~iH{t — to)] I which formally solves 
the Schrodinger equation. It follows immediately that 
a stationary state, characterized by H\(j)i) = Ei\(j)i), 
with a real energy eigenvalue Ei, performs a simple har- 
monic motion, i.e., \ipi{t)) = exp[— i£'i(t — to)]\'ipi{to)) = 
e:icp[—iEi{t — to)]\(f)i). Henceforth, we assume a denumer- 
able set of such eigenstates of the Hamilton operator. 

Having recognized already the Hamiltonian character 
of the underlying equation(s) of motion, the harmonic 
motion suggests to introduce what may be called oscil- 
lator representation for such states. More generally, we 
consider the expansion of any state vector with respect 
to a complete orthonormal basis, {1$^)}: 



(11) 



i 



where the generally time dependent expansion coeffi- 
cients are explicitly written in terms of real and imag- 
inary parts, Xi,Pi. Employing this expansion, allows to 
evaluate more explicitly the Hamiltonian function intro- 
duced in Eq. ©, i.e., U = {^H\^): 



n 



2 ^ 



['t>,\H\'P,){X,~iP,){X,+iPj) 



H{X„P,) 



(12) 



Choosing especially the set of energy eigenstates, {|0i)}, 
as basis for the expansion, we obtain: 



nix^.,p.)^Y.Y^P" + ^") 



(13) 



hence the name oscillator representation. The simple rea- 
soning leading to this result clearly indicates that {Xi, Pi) 
may play the role of canonical coordinates in the descrip- 
tion of a quantum mechanical object and its evolution 
with respect to the state space. - However, several points 
need to be clarified, in order to validate this interpreta- 
tion. 
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First of all, with {Xi,Pi) as canonical coordinates 
and H as Hamiltonian function, we verify that the 
Schrodinger equation is recovered by evaluating = 
\^i){Xi + iPi)/\/2 according to Hamilton's equations 
of motion: 



X, = 



dp 



It remains to demonstrate the compatibility of the no- 
tion of observable introduced here - as in classical me- 
chanics, cf. the discussion leading to Eq. ([2]) and there- 
after - with the one adopted in quantum mechanics. This 
concerns, in particular, the implementation of canonical 
transformations and the role of observables as their gen- 
erators. 



J2 {h^AXj + iPj) - {X, - iP,)hA , (14) 



3. Canonical transformations and quantum observables 



dnx,,p,) 

dX,, 



(15) 



where Hij 



Inserting these terms 

and using Eq. (jlip leads to \'i>) = —iH\'i'), as expected. 
Using "H in the special form given by Eq. (IT51) . we see 
that a zero mode with Eii = automatically leads to 
{Xir,Pi>) — constant. 

Secondly, the constraint C := = 1, cf. Eq. ([5]), 

becomes: 



C{X„P,) 



(16) 



Thus, the vector with components given by the canoni- 
cal coordinates {Xi,Pi), i = 1, . . . , A^, is constrained to 
the surface of a 2A^-dimensional sphere with radius \/2. 
This constraint obviously presents a major difference to 
classical Hamiltonian mechanics. 

Following the previous discussion in Subsection 2.1., 
it is natural to introduce also here a Poisson bracket for 
any two observables on the spherically compactified state 
space, i.e. real- valued regular functions F, G of the coor- 
dinates {Xi,Pi): 



{F.G} = ^ 



fdF dG 
^ dxi 



dF dG 



(17) 



cf. Eq. ([T]). - Then, as before, the Hamiltonian acts as 
the generator of time evolution of any observable O, i.e.: 



(18) 



In particular, it is straightforward to verify with the help 
of Eqs. (Il4])-(I15]) that the constraint, Eq. ([16]), is con- 
served under the Hamiltonian flow: 



(19) 



Therefore, it is sufhcient to impose this constraint, which 
implements the normalization of the quantum mechanical 
state, on the initial condition of time evolution. 



The Hamiltonian function has been introduced as ob- 
servable in the Eq. ()12p which provides a direct relation 
to the corresponding quantum observable, namely the 
expectation value of the self-adjoint Hamilton operator. 
This is an indication of the general structure to be dis- 
cussed now. 

Refering to Section HI. of Heslot's work [12] for details 
of the derivations, we summarize here the main points, 
which will be useful in the following: 
• A) Compatibility of unitary transformations and Pois- 
son structure. - The canonical transformations discussed 
in Section 2.1. represent automorphisms of the classi- 
cal state space which are compatible with the Poisson 
brackets. In quantum mechanics automorphisms of the 
Hilbert space are implemented by unitary transforma- 
tions, I*') = t/|*), with UU'' = U'<U = 1. This implies a 
transformation of the canonical coordinates here, i.e., of 
the expansion coefficients {Xi, Pi) introduced in Eq. ()11|) : 



I* 



X,- -t- iPj 



X' + iP' 



V2 



(20) 



Splitting the matrix elements Ui- 



($,|{7|$j) into 



real and imaginary parts and separating Eq. (j20p accord- 
ingly, using orthonormality of the basis, yields the trans- 
formed coordinates in terms of the original ones. Then, 
a simple calculation, employing the Poisson bracket 
defined in Eq. ([TTl), shows that {X[,P'^ = (5„ and 
{X'i,X'j} = {i^',P;} = 0, as before. The funda- 
mental Poisson brackets remain invariant under unitary 
transformations. More generally, this implies 27] that 
U{{F,G}) = {U{F),UiG)}, cf. Eq. (0). Thus,' unitary 
transformations on Hilbert space are canonical transfor- 
mations on the {X, P) state space. 

• B) Self-adjoint operators as observables. - Any in- 
finitesimal unitary transformation U can be generated 
by a self-adjoint operator G, such that: 



U = 1- iGSa 



(21) 



which will lead to the quantum mechanical relation be- 
tween an observable and a self-adjoint operator, replacing 
the classical construction in Section 2.1. In fact, straight- 
forward calculation along the lines of A) , splitting matrix 
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elements of G, with G*^ = Gij, into real and imaginary 
parts, shows that in the present case we have: 

X.^XH A, + 

Due to the phase arbitrariness - U and U ■ exp{i9), with 
constant phase 9, are physically equivalent - the operator 
G is defined up to an additive constant. This constant 
is naturally chosen such that the relation between an ob- 
servable G, defined in analogy to Section 2.1., and a self- 
adjoint operator G can be inferred from Eqs. 

G{X,,P,)^{^\G\^) , (24) 

by comparison with the classical result, Eqs. ©-Q- In 
conclusion, a real-valued regular function G of the state 
is an observable, if and only if there exists a self-adjoint 
operator G such that Eq. \24^ holds. - Note that all quan- 
tum observables are quadratic forms in the XiS and Pi's. 
This explains that there are much fewer of them than in 
the corresponding classical case. 

• C) Commutators as Poisson brackets. - The relation 
(1231) between observables and self-adjoint operators is lin- 
ear and admits 1 as unit operator, since = 1. 
Therefore, addition of observables and multiplication by 
a scalar of observables are well-defined and translate into 
the corresponding expressions for the operators. One 
may then consider the Poisson bracket (IT7)) of two ob- 
servables and demonstrate the important result [T2j : 

{F,G}^{^\-[F,Gm , (25) 
I 

with both sides of the equality considered as functions of 
the variables Xi, Pi, of course, and with the commutator 
defined as usual, [F, G] := FG — GF. This shows that 
the commutator is a Poisson bracket with respect to the 
(X, P) state space and relates the algebra of observables, 
in the sense of the classical construction of Section 2.1., 
to the algebra of self-adjoint operators in quantum me- 
chanics. 

• D) Normalization, phase arbitrariness, and admissible 
observables. - Coming back to the normalization con- 
dition ('5 1 'J) = 1, which compactifies the state space, 
cf. the constraint Eq. (jl6p . it must be preserved un- 
der infinitesimal canonical transformations, since it be- 
longs to the structural characteristics of the state space. 
By Eqs. an infinitesimal canonical transforma- 
tion generated by an observable G leads to C(Xi, Pi) 
C(X,',P/), with: 

/ fIG BG \ 

C{X[,P^ ^ CiX,,P,) + J2[—X,-—P,)5a 

j 3 3 

^0{5a^) . (26) 



Therefore, a necessary condition which observables must 
fulfil is the vanishing of the term cx 5a here, i.e. the 
invariance of the constraint under such transformations, 

C(X', P/) = C{Xi, Pi). It is not sufficient, since the prod- 
uct G1G2 of two observables - which fulfil this condition 
individually and, therefore, their product as well - does 
not necessarily represent an observable: the correspond- 
ing self-adjoint operators do not necessarily commute, 
i.e., generally we have (G'iG2)^ = G2G1 7^ G1G2. - In- 
cidentally, the condition of the vanishing second term on 
the right-hand side of Eq. (pS)) follows also more gener- 
ally, via Eq. (|24p. from the requirement that any observ- 
able G is invariant under an infinitesimal phase trans- 
formation Y^) — > 1^*) • exY>{i56), with constant 58, 

G{X[,P[) = G{Xi,Pi). Conversely, assuming this phase 
invariance of the observables, we recover that Hilbert 
space vectors differing by an arbitrary constant phase are 
indistinguishable and represent the same physical state. 

We note that any observable G with an expansion as 
in Eq. (1121) automatically satisfies the invariance require- 
ments of item D) above, the vanishing of the second term 
on the right-hand side of Eq. ([^5]). in particular. Explicit 
calculation shows: 

iC.O)-^{§-X,-§.p,)... ,27, 

3 

assuming that: 

G{P,,Xi) := (vE-IGI*) = i^G.,(X,-zPO(X,+iP,), 

i,3 

(28) 

and where G^ :— ($i|G|<i>j) — G*i, for a self-adjoint 
operator G. 

In conclusion, quantum mechanics shares with classi- 
cal mechanics an even dimensional state space, a Poisson 
structure, and a related algebra of observables. Yet it 
differs essentially by a restricted set of observables and 
the requirements of phase invariance and normalization, 
which compactify the underlying Hilbert space to the 
complex projective space formed by its rays. 

III. HYBRID PHASE SPACE ENSEMBLES 

So far, we have described the Hamiltonian formalism 
of classical mechanics and its generalization which cov- 
ers quantum mechanics, by adding more structure to the 
relevant state space. With the Hamiltonian equations 
of motion at hand, we could proceed to study the evolu- 
tion and direct coupling of classical and quantum objects. 
However, it is convenient to study the evolution of ensem- 
bles over the state (or phase) space instead liOj, JJjj . Last 
not least, this will allow to include quantum mechanical 
mixed states, thus generalizing beyond the above tacitly 
assumed pure states. 

In this section, we still neglect interactions between 
classical and quantum sectors of a combined system, the 
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study of which will lead us to truly quantum-classical hy- 
brids only in the next Section 4. - From now on, we will 
refer to the classical and quantum sectors as "CL" and 
"QM" , respectively. 

We describe a quantum- classical hybrid ensemble by a 
real-valued, positive semi-definite, normalized, and pos- 
sibly time dependent regular function, the probability dis- 
tribution p, on the Cartesian product state space canoni- 
cally coordinated by 2(n -I- A^)-tuples {xk,Pk', Xi, Pi); we 
reserve variables Xk,Pk, k = 1, . . . ,n for the CL sector 
(cf. Section 2.1.) and variables Xi,Pi, i = 1, . . . , N for 
the QM sector (cf. Section 2.2.). One or the other sector 
of the state space may eventually be infinite dimensional. 
A physical realization of such an ensemble can be imag- 
ined as a collection of representatives of the combined 
system with different initial conditions. 

In order to qualify as observable, the distribution p 
additionally has to obey the constraint induced by the 
extra structure of the QM sector of state space, see Sub- 
section 2.2.3. D). Evaluating the expectation of the cor- 
responding self-adjoint, positive semi-definite, trace nor- 
malized density operator p in a generic state j^*), Eq. dill) , 
we have p{xk,Pk; Xi, P^) := {^\p{xk,Pk)\^) , and: 

p{xk,Pk] Xi,Pi) = - '^pij{xk,Pk){Xi - iPi){Xj + iPj), 

(29) 

with pij{xk,Pk) ■= {^i\pixk,Pk)\^j) = p*i{xk,Pk)- This 
assures that p (or the marginal QM distribution obtained 
by integrating over the CL variables), as generator of a 
canonical transformation, does not violate the normal- 
ization constraint and phase invariance; in particular, it 
follows that {C,p} = 0, cf. Eq. (^7)) . 

Furthermore, positive semi-definiteness of p imposes 
constraints on any other observable G (g) of the QM 
(CL) sector, which can generate a canonical transforma- 
tion. Considering infinitesimal transformations in both 
sectors, cf. Eqs. ©-g]) and Eqs. ([22]) -(gH), we obtain 
p{xk,Pk;Xi,P,) p{x'^,p'^;P.-,Xl), with: 

p{xk,Pk;Xi,Pi) = p{xk,Pk\Pi,Xi) 

+{dx^p dp^g - dp,p dx^g)SacL 
+{dx,p dp.G - dp.p dxiG)SaQM 
+0(focL,(5aQM,^acL^Q;QM) -(30) 

Now, if and where the distribution p vanishes, also the 
first order terms cx ^acL and cx daQM have to vanish, 
since otherwise p can be made to decrease below zero by 
suitably choosing the signs of these independent infinites- 
imal parameters. 

This will be particularly relevant for the time evolution 
generated by a quantum-classical hybrid Hamiltonian, to 
be discussed in Section 4. 



A. The probability density and marginal 
distributions 

Finally, we remark that the relation between an observ- 
able in {X, P)-space and a self-adjoint operator, Eq. (|24l) . 
can be written as: G{Xi,Pi) = Tr(|*)(*|G), which 
shows explicitly the role of a QM pure state as one- 
dimensional projector, in this context. In order to illu- 
minate the meaning of the probability density p, we may 
then use the representation of p in terms of its eigen- 
states, p = J2j ^iid obtain: 

p{xk,Pk;X,,P,) = ^«;,(xfe,pfc)Tr(|*)(^'|j)(j|) 
i 

= Y^w.ixk^PkWm^ , (31) 

with <Wj <1 and J^j J ^i{<^xidpi)ujj{xk,Pk) = 1- 

We see that p{xk,PkT Xi, Pi), when properly normal- 
ized, is the probability density to find in the hybrid en- 
semble the QM state parametrized by Xi, Pi through 
Eq. PT|) . together with the CL state described by the co- 
ordinates {xk,Pk) of a point in CL phase space. 

The probability density allows to evaluate expectations 
of QM, CL, or hybrid observables in the usual way. Par- 
ticularly useful are also the marginal (or reduced) distri- 
butions defined by: 

Pci.{xk,Pk) ■■= 

/ YVjidXjdPj) p{xk,Pk;X„Pi) , (32) 

PQM{Xi,Pi) := J ni{dxidpi)p{xk,Pk;Xi,P,),{33) 

where Tn denotes a normalization factor, to be deter- 
mined shortly, and the integration in Eq. (15^ extends 
over the surface of a 2A^-dimensional sphere of radius V2, 
in accordance with Eq. (|16l) ; the integration in Eq. ([55)) 
extends over all the state space of the CL subsystem. 
The convergence of the integrals is assured by the pos- 
itive semi-definiteness and normalizability of p; the un- 
derlying assumption is that the CL subsystem occupies 
essentially only a finite region of its phase space, while 
the QM subsystem is constrained by the normalization 
of its state. 

More explicitly, using the representation given in 
Eg. (PT|) . we calculate for a state vector I^P), expanded 
according to Eq. (ITT|) : 

PCL{xk,Pk) ^^Wj{xk,Pk) ^{j\'^ii){'^i2\j) ■ hti2 , 
j ilM 

(34) 

with a remaining surface integral defined by: 
lab ■■= I n,(dA,dF,) 

J5S2n(V2) 

x{Xa+iPa){Xb-iPb) , (35) 
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and evaluated as follows: 



lab = (Ja^r^i /"n,(dX,dP,)(5(2-^(X,2 + P,2) 



Sab 



2JV-1 



s(r + V2){R~V2) 



Sab , 



(36) 



making use of isotropy and, in particular, replacing 
^a' + Pa by EaW + P^) I ^ = 2/iV Under the in- 
tegral; in the end, we employ 2iV-dimensional spherical 
coordinates, where VL2N denotes the spherical angle, and 
choose the normalization factor appropriately: 



N 



N 



(37) 



Thus, we obtain from Eq. ([34]) the expected simple result: 

Pci.{xk,Pk) ^^Wj{xk,Pk) , (38) 
j 

using completeness and orthonormality of the bases. 



B. Quantum-classical Poisson bracket and 
separability 

The result of the calculation in Eq. ([50)) suggests to in- 
troduce a generalized Poisson bracket, when considering 
observables defined on the Cartesian product state space 
of CL and QM sectors as follows: 



{A,B}, 



:= {A, B}cL + {A, B}qm 
dAdB dAdB\ 



Yi- 



E 



^dA dB 



dA dB 



(39) 



(40) 



for any two observables A, B. It is bilinear and antisym- 
metric, leads to a derivation-like product formula and 
obeys the Jacobi identity, since the right-hand side of 
Eq. (j40l) can be written in standard form as a single sum, 
after relabeling the canonical coordinates. 

Let us say an observable "belongs" to the CL (QM) 
sector, if it is constant with respect to the canonical coor- 
dinates of the QM ( CL ) sector. - Then, the generalized 
Poisson bracket has the additional important properties: 

• It reduces to the Poisson brackets introduced in 
Eqs. ([Ij and ([T7|) . respectively, for pairs of observ- 
ables that belong either to the CL or the QM sector. 



• It reduces to the appropriate one of the former 
brackets, if one of the observables belongs only to 
either one of the two sectors. 

• It reflects the separability of CL and QM sectors, 
since {A, i?}x = 0, if A and B belong to different 
sectors. 

The physical relevance of separability can be expressed 
as the following requirement: // a canonical tranforma- 
tion is performed on the QM ( CL ) sector only, then all 
observables that belong to the CL ( QM) sector should re- 
main unaffected. This is indeed the case, as we shall 
demonstrate directly by examining the behaviour of the 
reduced CL (QM) probability distribution under such 
transformations . 

Performing in the QM sector, for example, an infinites- 
imal canonical transformation on the integral of Eq. 
we obtain: 

Pci.{xk,Pk) ^ 



J, , n.idx'^dp;) p{xk,Pk;XiPl) - 

p[Xk,Pk;Xi H da. Pi — da) 



dP 



dX, 



= / p{xk,Pk;Xi,Pi) 



j {dx^p dp^G - dp^p dx,G)Sa + 0{Sc 



= PcLixk,Pk) + 0{Sa^) 
where we abbreviated P^^ 



(41) 



n,(dx,dP,) ^ / 

and used the well known invariance of the phase space 
volume element and of the constraint surface, by Eq. ()27p , 
together with Eqs. ([^^ ~ (|24l) : furthermore, the last equal- 
ity follows from the fact that the integral of a Poisson 
bracket of observables over QM state space vanishes: 



{A,B}qm = j {^\^[A,Bm 

Tr(|*)(*|i[i,B]) =Tr(i[i,B]) = , (42) 
i i 



using Eq. ()25l) , followed by a calculation similar to the 
one leading from Eq. ^ to Eq. ([5H|), via Eqs. ([32|) and 
Eqs. (IM1) - ([57|) . - In the present case, incidentally, we 
have that {A,B}qm = {p, G}qm = {p, G}x, since G 
belongs to the QM sector. 

Thus, we find invariance of pcL under infinitesimal 
and, hence, finite canonical transformations in the QM 
sector. Consequently, the expectation of any CL observ- 
able gcL, defined by: 

(ffCL) := rniidxidpi) gcL{xk,Pk)pcL{xk,Pk) , (43) 
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is invariant. - Similarly, one shows that pqm is invariant 
under canonical transformations in the CL sector and, 
thus, expectations of QM observables as well. 

Separability, as demonstrated here, has been a cru- 
cial issue in discussions of earlier attempts to formulate 
a consistent quantum-classical hybrid dynamics, see, for 
example, Refs. 0, 0, [H and references therein. 



C. Time evolution of noninteracting 
quantum-classical ensembles 

In order to illustrate another aspect of the separability 
of CL and CM sectors, as long as there is no interac- 
tion between them, we consider the time evolution of the 
probability distribution generated by the total Hamilto- 
nian function "Hs: 

H^ixk,Pk;X,,P,) :=HcL(a;fc,Pfc)+HQM(^^,-P^) , (44) 

where Hch denotes an assumed Hamiltonian function for 
the CL sector, while the Hamiltonian function Hqm for 
the QM sector has been detailed above, cf. Eqs. (|12p - 

m- 

Based on Hamilton's equations for both sectors 
and equipped with the generalized Poisson bracket of 
Eq. pop , we can invoke Liouville's theorem to obtain the 
evolution equation: 

-dtp = {p,ns}x ■ (45) 

Clearly, this equation admits a factorizable solution 
p{xk,Pk;Xi,Pi;t) = p{xk,Pk\t)- p{Xi, Pi;t), provided the 
initial condition has this property. No spurious correla- 
tions are produced by the evolution, which corresponds 
to {'Hcl,'Hqm}x — 0, by construction. 

In other words, the CL and QM sectors evolve inde- 
pendently, as if the respective other sector was absent, 
and maintain their classical or quantum nature, as long 
as they do not interact. 



IV. QUANTUM-CLASSICAL HYBRID 
DYNAMICS 

Following the preparations in Section 3., which con- 
cerned quantum-classical composite systems, however, 
without interaction of the CL and QM sectors, we pro- 
pose here the generalization to include such a hypothet- 
ical coupling and will study the consistency and conse- 
quences of such truly hybrid systems. 

This discussion will concern hybrid ensembles, or spe- 
cific hybrid states, and their dynamics. Given the gen- 
eralized Poisson bracket, introduced in Eq. ([5^ . we have 
to incorporate a hybrid interaction term T in the total 
Hamiltonian function "Hs, which will serve as the gener- 
ator of time evolution, as before. Therefore, we replace 
the definition of Eq. PH) . with "Hs = 'Hs(xfc,pfc; X^, P^), 



by: 

"Hs := Uchixk^Pk) +'HQM{Xi,Pi) +I{xk,Pk\Xi,Pi) . 

(46) 

For "Hs to be an observable, it is necessary that the hybrid 
interaction qualifies as observable, in particular. Further 
properties of "He will be detailed in due course. 

Then, the evolution equation is of the same form 
as Eq. (1451) : —dtp = {/3,'Hs}x , where, henceforth, 
the Hamiltonian function includes the interaction term 
I, unless stated otherwise. - Here, the positive semi- 
definiteness of p holds for the same reason as for the case 
of the classical Liouville equation, namely that the un- 
derlying dynamics is described by a Hamiltonian flow. 

A. Energy conservation 

Having proposed 'Hs as the generator of time evolu- 
tion, it also provides the natural candidate for the con- 
served energy of the hybrid system. Since Hj: is assumed 
not to be explicitly time dependent, we find, cf. with the 
general structure of Eq. (fT5)) : 

^ = {He,Hs}x -0 , (47) 

an immediate consequence of the antisymmetry of the 
generalized Poisson bracket. - Note that in the absence 
of a CL subsystem, this result reduces to the conserva- 
tion of the expectation of the QM Hamilton operator, 
as it should. More generally, in the absence of QM-CL 
interactions, the classical and quantum mechanical ener- 
gies simply add. 

B. Generalized Ehrenfest relations for hybrids 

Here we show that the Poisson structure built into 
the present theory of hybrid systems, in particular in 
the form of underlying Hamiltonian equations of motion, 
translates into generalizations of Ehrenfest relations for 
coordinate and momentum observables. 

We consider hybrid systems described by a generic 
classical Hamiltonian function and a quantum mechani- 
cal Hamiltonian operator, respectively: 

k 

Hqm ■■= ^ + V{X) , (49) 

where v{xi) = v{xi, . . . , Xn) and V denote relevant po- 
tentials, together with a self-adjoint hybrid interaction 
operator I{xk,Pk', X , P); note that symmetrical (Weyl) 
ordering is necessary, concerning the noncommuting op- 
erators X and P. We set all masses equal to one here, 
for simplicity, but will introduce them explicitly in the 
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particular case of coupled oscillators below. By Eq. (l24l) , 
this gives rise to the following Hamiltonian function H-^: 

fc 

+I{xk,Pk;X,,P,) , (50) 

when evaluated in a pure state j^*), invoking the oscillator 
representation of Eq. (Ilip . Correspondingly, we define 
coordinate and momentum observables, in the sense of 
our earlier construction in Section 2.2., pertaining to the 
QM subsystem: 

X{X„E,):^{^\X\^) , P(X„P,) (*|P|*) . (51) 

With these definitions in place, we proceed to determine 
the equations of motion by following the rules of Hamil- 
tonian dynamics. 

The equations of motion for the CL observables Xk , Pk 
are: 

Xk = {xk.HT.jx =Pk + dp^I{xk,Pk;Xi,Pi) , (52) 

Pk = {Pk^U^jx 

= -dx^,v{xi) - dx^T[xk,Pk]Xi,Pi) . (53) 

Similarly, we obtain for the QM variables Xi, Pi, which 
are not observables: 

Xi = {Xi,'HT,}x 

= dp,nQMiX,,P,) + dp^Iixk,Pk;Xj,P,) (54) 
= E,P,+dp^I{xk,Pk-Xj,Pj) , (55) 

h - {p^,n^}x 

= -5x.Hqm(^j,P,) - dx.I{xk,Pk;Xj,Pj){56) 
= -E,X,-dxMxk,Pk;X,,Pj) , (57) 

where Eqs. ([55)1 and (|57p follow, if the oscillator expan- 
sion is performed with respect to the stationary states of 
Hqm, cf. Eqs. (dUl-dni) in Subsection 2.2.2. - Notably, 
the Eqs. ([521), ^ together with Eqs. or to- 

gether with Eqs. ([55|) . (|57)) . form a closed set of 2(n + N) 
equations, where n denotes the number of CL degrees of 
freedom and N the dimension of the QM Hilbert space 
(assumed denumerable, if not finite). 

However, in distinction, the generalized Ehrenfest rela- 
tions for the QM observables X,P, defined in Eqs. (1511) . 



are obtained as follows: 
X = {X,Hs}x ={^,He}qm 

= -i{^S\[X,HQM + im 

= p-tmxj{xk,Pk;x,p)m , (58) 

p = {Pn^}x = {phs}qm 
= -?:(*|[P,i?QM + /]|*) 

= ^{^\r{x)\^)-tmpjixk,Pk;x,p)m,{59) 

where we used Eq. ((25|) , in order to replace Poisson 
brackets by commutators and the explicit form of Hqm, 
Eq. (P5|) : V' denotes the appropriate first derivative of 
the potential function V. - The Eqs. (l55])-(in5|) together 
with Eqs. (|52p - (j53p do not form a closed set of equations, 
since the expectation of a function of observables gener- 
ally does not equal the function of the expectations of 
the observables, as in Ehrenfest's theorem in quantum 
mechanics. [29| 

C. Bilinearly coupled oscillators 

We consider here a set of CL oscillators coupled bilin- 
early to one QM oscillator, choosing, for example: 

^CL := E i^kP' + '^-') ' (60) 

k 

Hqm ■= j^P"^ + ^^^X"^ , (61) 

/ := l^AfcXfc , (62) 
fc 

where we introduced masses mk,M, frequencies a;fc,f2, 
and coupling constants Afc. 

In this case, the equations of motion for the CL observ- 
ables together with the generalized Ehrenfest relations of 
the previous subsection reduce to a simple closed set of 
equations: 

= ' (63) 

Pk = -nikuji^xk - XkX , (64) 
X = jjP , (65) 

P = -Mfl^X -J2>^kXk , (66) 

fc 

with the QM observables X :== (*|^|*> and P := 
(\I'|P|vE'), cf. Eqs. ([51]). Here, the backreaction of QM 
on CL subsystem appears, as if the CL subsystem was 
coupled to another CL oscillator. 

In view of Eqs. we find that our theory 

passes the "Peres- Terno benchmark" test for interact- 
ing QM-CL hybrid systems which, so far, has been 
achieved only by the configuration ensemble theory of 
Hall and Reginatto H-Q. 
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DISCUSSION 



a driven harmonic oscillator. This yields: 



The proposed theory describing QM-CL hybrid sys- 
tems certainly raises a number of questions, some of 
which we address in the following. 



A. Classical-environment induced decoherence 

Well known studies of environment induced decoher- 
ence describe the effects that an environment of QM 
degrees of freedom has on the coherence properties of 
a QM subsystem coupled to it [3l|, [s^. In particular, 
the Feynman- Vernon or Caldeira-Leggett models and, 
more generally, models of quantum Brownian motion 
have been studied in this respect [111. Here we suggest 
to consider the situation where the QM environment is 
replaced by a classical one. We shall find that a CL envi- 
ronment similarly can produce decoherence in a generic 
model. 

For simplicity, we consider a QM object characterized 
by a two-dimensional Hilbert space, a "q-bit" , which is 
coupled bilinearly to a set of CL oscillators. The oscilla- 
tors are described by Eq. as before. The Hamilto- 
nian function of the q-bit presents the simplest example 
of the oscillator expansion of a QM Hamiltonian: 



(67) 



when expanding with respect to the energy eigenstates, 
{|</'i),l02)}, cf. Eqs. The bilinear QM-CL cou- 

pling is defined by: 



/ := 



(68) 



similarly as before in Eq. (|62|) : here E denotes an observ- 
able of the q-bit. 

In this case, the closed set of dynamical equations of 
motion becomes: 



(69) 
(70) 
(71) 

Pi = -E^Xi > AfeXfe , (72) 

Ox, ^ 



analogous to Eqs. ([5^ . ([55)) . ([55)) . ([57)) . and where we 
have: 

(vflSlvl/) = 1 ^ (0,|E|0,)(X,-zPO(^,+*^',) • (73) 

The Eqs. (IMt - lfTO]) are solved by employing the re- 
tarded Green's function for the equation of motion of 



(74) 

with the harmonic term x"^^ it) := afccos(ajfet) + 
hk sin(a;fei) and where the real coefficients ak and bk are 
determined by initial conditions. 

Furthermore, the Eqs. ([7T|) - ([72)) can be combined into 
second order form: 

Xi + E,j^Xi 

= -EU{t) E (Re(S,,)^, -MSy)^j/^j)), (75) 

introducing the matrix elements := ((/)i|S|^j) = E*j, 
the real and imaginary parts of which enter. Thus, 
we obtain a system of iV = 2 coupled oscillator equa- 
tions, where the coupling terms are nonlinear and non- 
Markovian through the function: 

m ■■= E^fe^'=(^) 

= eM"^(^) 



X (x,is)X,{s) + X,{s)X,{s)/{E,E, 



(76) 



by Eqs. (HB), GSj-dZlD. 

Let us reduce the above model to a crudely simplified 
version, neglecting presumably much of the rich dynam- 
ics described by Eqs. ([75l) -(|76 |) . - For sufficiently weak 
coupling, we drop the non-Mar kovian terms, i.e., terms 
oc . Furthermore, we choose En — E22 = and 
E12 — E21 = 1. This simplifies the equations to de- 
scribe two oscillators which are symmetrically coupled to 
each other through a periodic or, in the case of CL oscil- 
lators with incommensurate frequencies, quasi-periodic 
function ^. - Under the additional assumption of slow 
CL oscillators, i.e., with frequencies that are small com- 
pared to the ones of the QM oscillators, the resulting 
equations are solved by (i = 1,2): 



Xi = Ai cos{Vtit) + Bi sin(r2ii) 



(77) 



and Pi — Xi/Ei, the real coefficients Ai, Bi are deter- 
mined by initial conditions and, to leading non-vanishing 
order in ^, the characteristic frequencies are given by: 



Sli := El 



2{E,' - Ei 



, ^2 



Eo 



2{E,^~E,^) 



(78) 

Choosing, for illustration, initial conditions such that 
the expansion coefficients in |^) — X^i 2 l*^*) ("'^^ + 
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iPi)/y/2 (cf. Eq. pTjl ) are real at t = 0, we obtain the 
off-diagonal matrix elements of the corresponding density 
matrix p := in the form: 

{4>l\p\4>2) = {4>2\p\4>l)* = (^1 + tPl){X2 - ^P2)/2 

+0^) . (79) 

We note that there is a term cx i^^i contributing to the 
argument of each exponential, cf. Eqs. ([75)) . This can 
lead to decoherence by dephasing in the following way. 

If the nonnegative function {t) is sufficiently irregular 
(depending on the frequency distribution of environmen- 
tal oscillators), we may treat it as a random variable and 
average the result of Eq. ([79)) correspondingly. We con- 
sider the leading term, while the others can similarly be 
dealt with. Thus, writing VI2 - 0.i ^ 5E + (^I25E, with 
5E := E2 — El, we have to evaluate the dimensionless 
function /: 

f{t) := / dn P{n)e'^' , (80) 
Jo 

with fl = /25E, and where P represents the appropri- 
ately normalized distribution of the values of Jl. Now, 
there are continuous distributions, such that f(t) — 0, 
for t ^ oo] for example, a constant distribution over a 
finite range of £7, an exponential distribution, or a Gaus- 
sian distribution. Under these circumstances, the leading 
term (similarly the others) gives a decaying contribution, 
i.e. oc f{t), to the off-diagonal density matrix element 
{4>i\p\4>2), after averaging. 

This indicates a decoherence mechanism which is effec- 
tively quite similar to "fundamental energy decoherence" , 
which has been reviewed recently in Ref. |19i] . 



set of canonical variables, for example, (^i, -Pi)i=i,...,7v 
in the closed set of dynamical equations ([5^ - ([57)) . As 
long as no approximations are applied to these equa- 
tions, their solutions allow to evaluate exactly all quan- 
tities which reflect the fluctuations associated with a 
pure quantum state |^'), such as the typical variance 
AX^ := {^X'^\^) - This follows from the 

fact that all quantum observables can be expanded in the 
oscillator representation, recall Eqs. (|TT|) . (|24)) . (p8)) . with 
the expansion coefficients given by the solutions of the de- 
terministic equations. Thus, for example, AX^ becomes 
a function of the canonical variables. 

The QM variables do not fluctuate in a given pure 
state. By the QM-CL Poisson brackets and ensuing equa- 
tions of motion ("Hamilton's equations") they are cou- 
pled to the CL variables (observables) which, therefore, 
do not show fluctuations either. 

However, this admits the possibility that the initial 
conditions of the hybrid dynamical equations, in partic- 
ular for the QM subsystem, are determined by the fluc- 
tuating outcome of a certain preparation / measurement. 
In this case, if only statistical / conditional information 
is available about the initial state of the system, the 
classical observables, generally, will reflect corresponding 
fluctuations. For example, we can evaluate a correlation 
function of CL observables to find: 

{XaXb) ■■= J Ui{dxidpi) XaXtPCL{xk,Pk) ^ {Xa){xb) , 

with the help of the reduced distribution pcL introduced 
in Eq. p2p . This distribution function is determined by 
the solution of the Liouville equation for the full density p 
of the interacting hybrid system, cf. Eq. (|45|) : it could be, 
furthermore, conditioned by a selected outcome of some 
quantum measurement (s) specifying the initial state. 

C. Hybrid observables, separable interactions and 
QM-CL Poisson brackets 



B. Quantum-classical backreaction 

Quantum-classical backreaction, in particular the ef- 
fect of quantum fiuctuations on the classical subsystem, 
has always been an important topic for various propos- 
als of quantum-classical hybrid dynamics and its ap- 
plications. This concerns improvements of approxima- 
tion methods and applications, for example, in "semi- 
classical gravity" studying the effects of quantum fiuc- 
tuations of matter on the classical metric of spacetime; 
see Refs. 0, @, H, [11] with numerous references to related 
work. 

Our formalism consistently incorporates all quantum 
fluctuations, even if they are not explicitly visible, un- 
like in many approaches where fluctuations are added 
by hand, in some approximation. Presently, the quan- 
tum dynamics is treated exactly in terms of a complete 



It is a common feature of either QM or CL systems 
that particular forms of interaction among subsystems 
allow to separate degrees of freedom into noninteracting 
subsets. Generally, this is associated with the existence of 
symmetries of the compound system, such as translation 
or rotation invariance. 

A recent study of a translation invariant harmonic in- 
teraction between a QM and a CL particle revelled that 
- according to the hybrid theory proposed by Hall and 
Reginatto - there arises an irreducible coupling between 
center-of-mass and relative motion Q. This is contrary 
to what happens if both particles are treated as either 
classical or quantum mechanical and has been traced to 
the inherent nonlinearity of their proposal. The action 
functional, from which the equations of motion are de- 
rived, "knows" which variables belong to the QM and 
CL sectors, respectively, and mixing them by coordinate 
transformations produces the coupling. Consequently, 
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such a system has been proposed as a prospective testing 
ground, where their theory could be falsified experimen- 
tally 0. 

This issue can also be examined in the light of the 
present linear hybrid theory. - Specializing the system 
of bilinearly coupled oscillators of Subsection 4.3. as fol- 
lows: 

nci.:^^p\ Hqm:=^P^ I := X{x ■ 1 - X)^ , 

(81) 

we reconsider the example of Ref. [9] ; here 1 denotes the 
unit operator on the Hilbert space of the QM subsystem. 

As before, the Hamiltonian function which gener- 
ates time evolution of the composite system, ■— 
Hch + {'^\{Hqu + is conserved by construction, 

dUs/dt = 0, see Subsection 4.1. Furthermore, we know 
from Subsection 4.3. that the generalized Ehrenfest equa- 
tions for bilinearly coupled oscillators in terms of the CL 
observables, here x and p, and of the QM observables, 
X and P := {'^'\P\^), form a closed set, cf. 

Eqs. These equations of motion are nothing 

but Hamilton's equations for the "classical" Hamiltonian 
function: 

ni{x,p;X,P):=^p^ + ^P' + X{x-X)^ , (82) 

which implies that also is conserved, dH^/dt = 0. 
Then, it follows that the energy carried by quantum fiuc- 
tuations is separately conserved: 

^^{2Miip')-ip)')+H{x')-{xr)) 

= ^(Hs-H|l)=0 , (83) 

where (...) = (^f] . . . j^f). Note that all mixed QM-CL 
terms cancel and no energy is transferred between the 
(fluctuations of the) QM and the CL subsystem. 

We observe that the "classical" Hamiltonian is sep- 
arable by transforming the variables x,p and X.P to 
center-of-mass and relative variables: 

s := {MX + mx) /a , ps := P + p , (84) 
r:^X~x , p^:=f,(±P^l.p) , (85) 

with the total and reduced masses defined by cr := M+m 
and /i := M'm/{M + m), respectively. Thus, we find: 
= {p^ /2(t) + {p^/2^i) + Ar^, not surprisingly. 
At this point, it seems worth while to assess the char- 
acter of these transformations with respect to the fun- 
damental QM-CL Poisson brackets, defined in Eqs. p9p - 
(PU]) (Subsection 3.2.), on which our theory is based. We 

have {x,p}x = 1 and {X,P}x = {[X,P]/i) = 1, us- 
ing Eq. ([25]) : similarly, we find that the brackets of all 
other pairs of these variables vanish. Thus, we may con- 
sider x,p and X, P as two pairs of canonical phase space 
coordinates. Furthermore, one may check that also the 
two pairs of center-of-mass and relative variables, s^ps 



and r, Pr, respectively, form pairs of canonical coordi- 
nates under the QM-CL Poisson brackets. Therefore, the 
transformations (j84p - (|85|) are consistent canonical trans- 
formations. 

An immediate consequence is that the "classical 
Hamiltonians" describing center-of-mass and relative mo- 
tion are separately conserved as well. 

It must be emphasized that the separation of the "clas- 
sical" degrees of freedom, x,p and X,P, or of the corre- 
sponding center-of-mass and relative variables, from the 
full set of canonical variables x,p and Xi,Pi, is an acci- 
dent of the harmonic interaction. Independently of the 
hybrid coupling between a classical and a quantum me- 
chanical particle, cf. Eqs. (I5T|) . any other but constant, 
linear (in one dimension), or harmonic translation in- 
variant coupling would not allow such separation, even if 
both particles were treated quantum mechanically. 

In general, the "mean field" variables X = {X) and 
P = (P) will always couple to "correlation functions", 
such as (X^), {XP + PX), {P^), or more comphcated 
ones, depending on the kind of interaction. This phe- 
nomenon of quantum mechanics is not particular to hy- 
brid dynamics. 

While the separation of degrees of freedom in the case 
of translation or rotation invariant potentials in quantum 
mechanics can always be completed at the operator level, 
the hybrid dynamics presented here necessitates the con- 
sideration of "classical" canonical variables on which to 
perform any canonical transformations consistently with 
the underlying Poisson bracket structure. This seems to 
limit separability to certain potentials, as we have just 
seen. 

We conclude that a composite system of a QM and 
a CL particle with harmonic translation invariant inter- 
action, or some analogue of this, does not allow to ex- 
perimentally falsify our formulation of hybrid dynamics. 
We find no coupling between relative and center-of-mass 
motion, contrary to the proposal of Refs. However, 
anharmonic interactions need to be studied in this con- 
text and may lead to experimentally accessible signatures 
of the linear QM-CL hybrid dynamics. 



D. The classical x almost-classical algebra of 
hybrid observables 

In Subsections 2.1. and 2.2., we introduced the no- 
tions of classical and quantum observables, respectively, 
relevant for the considerations of this paper. 

Furthermore, in Subsection 3.2., we introduced the 
fundamental QM-CL hybrid Poisson bracket, {A, B} x := 
{A,B}cL + {A,B}qm- Following Eqs. (I39l)~(l40l), we 
pointed out three important properties of this bracket, 
last not least related to separability. However, we tacitly 
assumed that a fourth case would not need further men- 
tion, which can arise for two genuine hybrid observables: 

• Let A = A{xk,Pk;X,,Pi), B = B{xk,Pk; Xi, P,) 
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be both hybrid observables, i.e., both are quadratic 
forms in the X^'s and P^'s and both are not com- 
pletely independent of the x^s and pk 's. If, further- 
more, one observable, say A, depends on any pair 
of canonical variables, say x = xw and p = pk' , and 
B also depends on x or p, then the "classical part" 
of the bracket, {A, -B}cl, generates terms which do 
not qualify as observable with respect to the QM 
sector. 

Such terms are of the general form: 

M,^j^,,^j,{xk,Pk){X, - iPi){Xj + iPj) 

x{X,. -iP,,){X,, +iP,,) 

= 4 ^ {^^^){^^^')^'h^J^,,^J,{xk,Pk) 

x(cf>,|M.) , (86) 

where we used the oscillator expansion, Eq. ([TT|), and: 

M,,,,,,{xk,Pk) ■■^-JZ[-9^-d^ - d d 

using the related expansion for observables A and B, cf. 
Eq.dM])- 

Generally, iterations of such brackets will implicitly 
contribute to the solution, p = p[xk,Pk] Xi^ Pi)^ of the 
evolution equation, —dtp — {p, "Hsjx , in the presence 
of a true hybrid coupling, cf. Section 4. Thus, multi- 
ple factors involving the state vector j^*) and its adjoint, 
or multiple pairs like {Xi — iPi){Xj -\- iPj), will enter. In 
this way, evolution of hybrid observables, of the density p 
in particular, can induce a structural change: while con- 
tinuing to be CL observables, they do not remain QM 
observables (quadratic forms in X^'s and Pi's). They fall 
outside of the product algebra generated by the observ- 
ables to which we confined ourselves, so far. 

We note that the assumption of a product algebra cov- 
ering the observables of a hybrid system was essential for 
the no-go theorem put forth in Ref . [IJl , which ruled out 
a class of hybridization models. However, this assump- 
tion can be criticized as being too restrictive from the 
QM point of view Q . - Here we assume: 

• The algebra of hybrid observables is closed under 
the QM-CL Poisson bracket { , }x operation - a 
physical hypothesis. 

Refering to the phase space coordinates (Xi,Pi), we 
define an almost- classical observable as a real- valued reg- 
ular function of pairs of factors like [Xi — iPi){Xj + iPj), 
such as in the left-hand side of Eq. (1551) . subject to the 

constraint: C{X„P^) = + Pi^) = 1- This nor- 

malization constraint, cf. Eq. ([T6|) in Subsection 2.2.2., 
is preserved under the evolution, since {C,'Hs}x = 0, in 
the presence of QM-CL hybrid interaction as defined in 
Section 4. 



According to this definition, QM observables 
(quadratic forms in phase space coordinates, cf. 
Section 2.2.) form a subset of almost-classical ob- 
servables which, in turn, form a subset of classical 
observables (real- valued regular functions of phase space 
coordinates, cf. Section 2.1.). 

Furthermore, we may now say that members of the 
complete algebra of hybrid observables, generally, are 
classical with, respect to coordinates {xk,Pk) and almost- 
classical with respect to coordinates [Xi,Pi). 

This leads us to speculate about a physical conse- 
quence of the enlarged classical x almost-classical alge- 
bra for interacting QM-CL hybrids, as illustrated by the 
following Gedankenexperiment. 

Consider a quantum together with a classical object 
subject to a transient hybrid interaction. As long as the 
hybrid interaction is ineffective, both objects evolve in- 
dependently according to Schrodinger's and Hamilton's 
equations, respectively. However, once they form an in- 
teracting hybrid, the corresponding phase space density 
changes from a factorized form, in absence of any ini- 
tial correlation, to become an almost-classical/classical 
hybrid observable. Even if the hybrid interaction even- 
tually disappears, the density possibly maintains a mixed 
almost-classical/classical character. This agrees with the 
general structure of the evolution equation, yet needs to 
be understood in detailed examples. 

This outcome contradicts naive expectation that quan- 
tum and classical objects evolve separately in quan- 
tum and classical ways, after any hybrid interaction has 
ceased. - Two possibilities come to mind. Either persis- 
tence of the almost-classical/classical character is a phys- 
ical effect accompanying QM-CL hybrids, if they exist. 
Or our description needs to be augmented with a reduc- 
tion mechanism by which evolving observables return to 
standard QM or CL form (cf. Section 2.), following a 
hybrid interaction. Both possibilities seem quite inter- 
esting in their own right. We reserve this topic for future 
study. [Hi 



E. Hybrid dynamics and Wigner function approach 

Suppose a physicist unfamiliar with quantum mechan- 
ics were presented with the general equations of motion, 
Eqs. ((5^ - ([57)) (plus normalization constraint, Eq. (HH])). 

- We know that these equations present independent CL 
and QM sectors, in the absence of a hybrid interaction. 

- However, he/she would naturally interpret them to 
describe the dynamics of a composite CL object, with 
part of its phase space compactified, due to the con- 
straint. Thus, he/she finds perfectly local dynamics. In 
fact, our knowledge of nonlocal features can be traced 
to the definition of the canonical coordinates and mo- 
menta Xi,Pi, introduced by the oscillator representa- 
tion, Eq. (HI]), since: XjV2 = Re / dg <i>*{q)^{q) and 
Pj/A/2 = ImJ dg $*(q)\I'(g). Therefore, spatially nonlo- 
cal (and probabilistic) features enter by reference to the 



14 



QM wave function. [36{ 

In view of this, it might be surprising that our pro- 
posed hybrid dynamics passes the set of consistency re- 
quirements, cf. Section 1., in particular the requirement 
of conservation and positivity of probabihty, as we have 
seen. 

This must be contrasted with the problems that arise if 
one maps the QM sector "locally" on a would-be classical 
phase space by using the Wigner function approach and 
the corresponding version of the von Neumann equation. 
The latter differs from the classical Liouville equation by 
a series of corrections in powers of h which, in turn, incor- 
porate nonlocal features. It is well known that they spoil 
the interpretation of the Wigner function as a genuine 
probability distribution on phase space, since it gener- 
ally does not remain positive semidefinite, see Ref. [s^ 
for a comprehensive review on probability issues. 

Conversely, the dynamics of classical phase space dis- 
tributions, typically described by the Liouville equation, 
can be presented in quantum mechanical and, in partic- 
ular, in path integral form [l3, [HI- Again, the result- 
ing would-be quantum mechanical density matrix corre- 
sponding to a classical probability distribution is, gener- 
ally, not positive semidefinite. 

In both cases, the problem is caused by the interme- 
diate Fourier transformation, which apparently is not 
suited to represent the nonlocality properties appropri- 
ately, when formally relating phase space to Hilbert space 
and vice versa. Therefore, the Wigner function approach, 
which one could be tempted to employ, in order to sys- 
tematically reduce part of a composite QM system to a 
CL subsystem, thus defining a QM-CL hybrid, unfortu- 
nately violates the positivity of probability requirement. 

Similar problems have been encountered in Ref. 
where the QM— ^CL reduction is attempted via coherent 
"minimum uncertainty" states. 

In distinction, the oscillator representation allows to 
circumvent this difficulty, at the expense of introducing 
the phase space structure in an abstract way. [39| 



dom are directly coupled to each other. An important 
guideline has been to satisfy the complete set of consis- 
tency conditions mentioned in Section 1. and fulfilled, so 
far, only by the configuration space ensemble theory of 
Hall and Reginatto [6-8], while all earlier attempts failed 
in one or the other point. However, our linear theory de- 
viates from their nonlinear theory in that no 'spurious' 
coupling between center-of-mass and relative motion is 
found for a two-body system with a harmonic transla- 
tion invariant potential. 

This latter issue, quantum-classical backreaction, 
classical-environment induced decoherence, and comple- 
tion of the algebra of hybrid observables have been dis- 
cussed in Section 5., while further interesting topics are 
left for future work. These include: the hypothetical role 
of hybrid dynamics in measurement processes, seen as the 
interaction between a classical apparatus and a quantum 
object according to the Copenhagen interpretation, and 
the effect of classical and quantum degrees of freedom on 
entangled and classically correlated states, respectively, 
through hybrid interactions. 

On the technical side, since hybrid dynamics in the 
present formulation leads to a Liouville equation, as we 
discussed, the superspace path inte gral we have recently 
devised can be readily adapted to it |lol.[Tl|. This may be 
particularly interesting for applications in which hybrid 
dynamics is considered as an approximation scheme for 
complex quantum systems. 

In a more speculative vein, one could wonder about the 
essential difference between quantum and classical state 
spaces seen here, respectively, in the presence and ab- 
sence of curvature (cf. Section 2.2.2.). Does a properly 
understood classical limit of quantum mechanics, which 
possibly helps with the measurement problem [40|, with 
a commented extensive list of references given in the last 
one., require a dynamical treatment of this curvature? 
Conversely, is the hypothetical emergence of quantum 
mechanics from deterministic dynamics related to a dy- 
namical structure of phase space? 



VI. CONCLUDING REMARKS 
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